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The Kondo effect involves the formation of a spin singlet by a magnetic impurity and conduction 
electrons. It is characterized by a low temperature scale, the Kondo temperature, Tk, and an 
associated long length scale, £k = Kvf / '(fcsTjf) where vf is the Fermi velocity. This Kondo length 
is often estimated theoretically to be in the range of .1 to 1 microns but such a long characteristic 
length scale has never been observed experimentally. In this review, I will examine how £r- appears as 
a crossover scale when one probes either the dependence of physical quantities an distance from the 
impurity or when the impurity is embedded in a finite size structure and discuss possible experiments 
■ that might finally observe this elusive length scale. 

d ■ I. INTRODUCTION 



00 



The Kondo model gives a simplified description of the interaction of a single magnetic impurity with the conduction 
electrons in a host, with Hamiltonian: 



H » N = E *l%<^ + J ^mp ' 4(r = °)- C 1 ) 
Here is the impurity spin operator and S e \ is the electron spin density at position 

% i S el (r) = &Jf* P {?) (2) 

where ^> a (r) annihilates a conduction electron of spin a at r and repeated spin indices are summed over. We will 
focus on the case S- im p = 1/2 although much of the discussion here carries over to higher spin. The dispersion relation 

£h , will often be assumed to be that of free electrons: 
O 
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although that is not crucial. (For a general review of the Kondo model, see Ref. [l|.) Essentially the same model 
^\ ' can be used to describe a quantum dot in the Coulomb blockade regime in a spin 1/2 state, interacting with one or 
more metallic leads (which can be formed by a 2 dimensional electron gas or by quantum wires in a semi-conductor 
inversion layer). 

The dimensionless parameter which measures the strength of the Kondo interaction is Jv = Xq where v is the 
t-H ■ density of states at the Fermi energy. For the free electron dispersion in D dimensions, this is: 

OS'. »D = k^/icDVp) (4) 

O ■ 

where C3 = 2-7T 2 , C2 = 2-7T and c\ — it. Typically Jv <C 1 suggesting that perturbation theory could be useful. However, 
. ! perturbation theory encouters infrared divergences at low energy scales. In general, an n th order term in perturbation 
theory for some physical quantity characterized by an energy E is proportional to [Jv \n(D/E)] n , where T> is an 
ultra-violet scale of order the bandwidth or Fermi energy. It is found that the leading logarithmic divergences can be 
summed by expressing perturbation theory in terms of the renormalized coupling constant at scale E: 

\{E) w A + Xl \n(V/E) + ... (5) 

Note that X(E) increases in magnitude, as E is lowered, assuming it is initially positive (antiferromagnetic) . X(E) 
obeys the renormalization group equation: 

' /A /3(A) = -X(Ef + {l/2)X{Ef + ... (6) 



d\nE 

where /3(A) is the ^-function. Keeping only the quadratic term in this RG equation, we find 



A ^ - l-A lnW) W 
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Note that Ao, the bare coupling, is the value of the renormalized coupling at the ultraviolet scale, T>. As the energy 
scale is lowered, A becomes of 0(1) at the Kondo energy scale: 

k B T K =Ve- 1/x ". (8) 

This renormalization group equation, is derived by integrating out high energy Fourier modes, reducing the 
effective bandwidth from O(D), in energy units, to O(E). Once the bandwidth becomes narrow (as happens at low 
energies) the width in wave- vector is related to the width in energy by -E = vp%\k — Uf\ where vf is the Fermi velocity. 
Thus we may equally well describe X(E) as the effective coupling at wave-vector scale E/(Hvf) as at energy E. From 
this perspective, k B Tx / \hvp) is the characteristic wave-vector scale at which the effective Kondo coupling becomes 
large. Thus it is natural to introduce the Kondo length scale: 

i K = hv F /(k B T K ) (9) 

(We henceforth set k B and h to one.) It seems appropriate to think of the effective Kondo coupling as growing at large 
distances, becoming large at length scales of order Equivalcntly, we might expect physical quantities depending 
on a length scale r to be scaling functions of the ratio r/£x rather than depending on r and £rt separately. 

There is an interesting analogy here with Quantum Chromodynamics (QCD) the theory of the strong interactions 
in high energy physics. The renormalization group equation for the effective QCD coupling, gQCD, (describing the 
interactions of gluons with themselves and with quarks) is the same as Eq. © at quadratic order. It also gets large 
as the energy scale is lowered, and small as it is raised. In this case, one switches back and forth from energy units 
to momentum units using the velocity of light, c. The characteristic energy scale, Aqcd where gQCD becomes 0(1) 
is of order 1 GeV, the mass of the proton, and the corresponding length scale is or order the Compton wavelength 
of the proton, 1 Fermi or 10 -15 m. If the quarks inside a proton are probed with high energy (short wavelength) 
photons they appear nearly free. On the other hand, they exhibit confinement (with an interaction growing linearly 
with separation) at long distances. In high energy physics it is commonplace to go back and forth freely between 
energy and length units using c (and K). However, there are well understood pitfalls in taking this picture of the 
effective coupling getting weak at high energies and strong at low energies, too literally. (These pitfalls occur both in 
the energy and distance picture, and do not seem to be particularly related to the distance viewpoint.) In this review 
article, we will be concerned with the validity of the corresponding picture for the Kondo model. 

The (5-function form of the Kondo interaction implies that only the s-wave harmonic interacts with the impurity 
so that the model becomes fundamentally 1-dimensional. After linearizing the dispersion relation about the Fermi 
surface, the low energy effective Hamiltonian becomes: 



H = ^l dr 
2tt 



VFXS imp ■ S el (0). (10) 



Here r is the radial coordinate and ipL/R represent incoming and outgoing waves, with the boundary condition: 

Mo) = Mo)- (ii) 

These are defined in terms of the s-wave part of the 3D electron annihilation operator, ty(r) by: 

*(r) = -±- [e' ikpr Mr) ~ e^ r ^ R {r)] +... (12) 

where the . . . represents higher spherical harmonics. Note that we have normalized the fermion fields as in Ref. [2| so 
that: 

{iPl(r),Mr')}=2ri(r-r'). (13) 

In the limit of zero Kondo interactions, Ao = 0, this describes a relativistic Dirac fermion with the Fermi velocity 
playing the role of the velocity of light. In such a model it is natural to go back and forth between energy and length 
units using vf ■ Note that it is crucial to this estimate of £k that only one channel (the s-wave) couples to the Kondo 
impurity, allowing a mapping into a ID model. While this happens in a variety of circumstances, including the case 
of a quantum dot coupled to 2DEG's, there are also important cases where it can fail, which will be discussed later. 

An intuitive picture of this Kondo length scale can be obtained from considering the low energy strong coupling 
behavior of the model. This is most easily understood from a 1 dimensional tight-binding version of the model with 
Hamiltonian: 

oo 

H = -*X>}^+1 + h - c -) + J ^imp ' (14) 

j=0 
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The strong coupling limit is easily understood. When J ^> t > 0, one electron gets trapped at the origin to form 
a singlet with the impurity spin. Note that this "Kondo screening" actually corresponds to the formation of an 
entangled state between the impurity spin and one conduction electron. The other electrons can do whatever they 
want except that they cannot enter or leave the origin since this would break up the Kondo singlet, costing an energy 
of O( J). They effectively feel an infinite repulsion from the origin, corresponding to a 7r/2 phase shift. This is simply 
a boundary condition on otherwise free electrons. 

While the strong coupling limit is trivial, we are actually interested in the case of weak bare coupling, at low 
energies and long distances, where the effective coupling becomes strong. It is known from various approximate 
and exact calculations that the low energy physics (at E <C Tk) is described by a local spin singlet state and non- 
interacting electrons with a 7r/2 phase shift. On the other hand, the physics at intermediate energy scales, of 0(Tk) 
is complicated. To form a spin singlet with an S=l/2 impurity, only one electron is needed and an intuitive picture 
is that one electron is "removed from the Fermi sea" for this purpose. However, unlike the simple case of large bare 
coupling, it would be quite wrong to think that this electron is localized at the origin. The natural length scale over 
which we may think of this electron's wave- function being non- negligible is £k- Such a naive picture must be used 
with caution. At best, it is valid only at long distances and low energies. If we probe the screened impurity with 
a long wavelength probe, this picture may apply. At shorter distances, it certainly breaks down. While this is only 
an intuitive picture, it nonetheless seems reasonable that £k will appear as a characteristic scale in any distance 
dependent physical property of a Kondo system. The nature of the crossover at ^k is the subject of this review. 

Nozieres' local Fermi liquid theory is well-known to provide a powerful way of studying the behaviour of distance- 
independent quantities at low T <C Tk- It turns out to also be useful for studying distance-dependent quantities at 
large distances r 3> £k- I n this approach, the screened impurity spin is eliminated from the effective Hamiltonian 
which, in lowest approximation, just consists of non- interacting electrons with a 7r/2 phase shift. Then, interactions 
at the origin are added for these phase shifted electrons, which originate from virtual excitations of the screened 
impurity complex. We may simply list all possible interactions at the origin allowed by symmetry. Assuming exact 
particle-hole symmetry, there is only one important leading irrelevant interaction: 



Here the current operator is defined as: 



"int = ~ ( 15 ) 



J L (r) = ^ L (r)^ L (r). (16) 



■^int con t ams a dimension 2 operator, so its coupling constant must have dimensions of inverse energy. On general 
grounds, we expect that this coupling constant will be of order the crossover scale, 1/Tr-. The factor of 1/6 in Eq. 
(TT51) is just a matter of convenience. With this normalization, a calculation of the impurity susceptibility, to first 
order in this interaction gives: 

X(T) - j^. (17) 

[This factor of 1/4 in x(T) is conventional since in the high T 3> Tk limit we obtain the free impurity result, 
x(T) — > 1/(4T).] We emphasize that various other definitions of Tk are in common use, differing by factors of 0(1). 

Perhaps surprisingly, distance dependent aspects of Kondo physics are much less well-studied than energy dependent 
ones. At a theoretical level this may be a consequence of the fact that it is very difficult to get distance dependence 
from two of the most powerful methods of studying the Kondo effect: the Numerical Renormalization Group (NRG) 
method, introduced by Wilson, 3 and the Bethe ansatz (BA) solution, discovered by Andrei^ and Weigmann^. The 
difficulty with NRG lies in the "logarithmic discretization" or the fact that an effective ID tight binding model 
is introduced in which the hopping parameter decays exponentially with distance from the impurity. This is, in 
fact, only an approximation to the full problem and one which seems to fail to keep track of length-dependence 
properly. However, more recently a way around this difficulty is being exploited^ which involves a more complicated 
and numerically costly version of NRG. The difficulty with BA seems to be that while exact wave- functions are 
calculated, they are given in such a complicated form that it is very challenging to calculate any Green's functions 
with them. Thus much of our understanding of length dependence comes from perturbation theory and analytic RG 
arguments, various mean field theories, exact diagonalization of short systems and Density Matrix Renormalization 
Group (DMRG) numerical results which are also restricted to fairly short systems (up to 32 sites). The reasons why 
this exponentially large scale, £k = v f/Tk, has never been seen experimentally are probably that the associated 
crossover effects at this length scale are rather weak and subtle and that the basic Kondo model may not be adequate 
to describe the experimental systems used. Some inprovements to the model that might be neccessary are: including 
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charge fluctuations (as in the Anderson model) taking into account a finite density of magnetic impurities, taking into 
account a finite density of non-magnetic impurities and taking into account electron-electron interactions even away 
from the impurity location. Thus searches for effects at this length scale represent both an experimental challenge to 
find sufficiently ideal systems and also an opportunity to study the limits of validity of the basic model. 

We will be concerned with two types of length dependence. The first type, analysed in Sec. UH involves a single 
impurity in an infinite host, described by the Hamiltonian of Eq. (JTJ. We consider various observables as a function 
of distance from the impurity. The first one f Sub-Section IPTA1 is the Knight shift, measurable in nuclear magnetic 
resonance (NMR) experiments. This is simply the magnetic polarization of the electrons as a function of distance 
from the impurity, in linear response to a magnetic field (applied to both the impurity and the conduction electrons). 
The second fSub-Section III B]) is the charge density (Friedel) oscillations, as a function of distance from the impurity. 
These could be observed by scanning tunneling microscopy (STM) for a magnetic impurity on a metallic surface. 
The third (Subsection III CP is the ground state equal time correlation function of the impurity spin and electron 
spin density at distance r from the impurity. While probably not experimentally observable, this has a rather direct 
interpretation in terms of measuring the spatial probability distribution for the electron forming the spin singlet with 
the impurity. The second type of length dependence, analysed in Sec. I1II| occurs in mesocopic samples containing a 
single localized S=l/2 where the size of some part of the device is comparable to £k- Here we consider four different 
situations. Two of them involve a S=l/2 quantum dot coupled to a finite length one-dimensional quantum wire. In 
the first case, (Subsection IIII A[) this wire is closed into a ring and we consider the persistent current through it as 
a function of the ring length and magnetic flux. In the second case, (Subsection IIII Bp the ring is straight with the 
quantum dot coupled to one end and the other end open. In Subsection IIII CI we discuss an S=l/2 quantum dot 
embedded in the centre of a finite length quantum wire which is tunnel-coupled to leads. In Subsection IIII Dl we 
contrast these strictly ID models with the case of a magnetic impurity inside a three dimensional magnetic sample 
containing non-magnetic disorder: the Kondo box models Given certain simplifying assumptions about the disorder, 
the characteristic Kondo length scale is much smaller than vf/Tk- It turns out that this is not just a consequence 
of the 3 dimensionality but also of the disorder. In an ideal 3D Kondo box £r- would again be the characteristic 
length, as we discuss. Electron-electron interactions away from the impurity are ignored in all sections except IIII BL 
a quantum dot at the end of a finite quantum wire, where they are included using Luttinger liquid techniques. Sec. 
IIVI contains conclusions. 



II. IMPURITY IN AN INFINITE HOST 



In this section, we consider the Hamiltonian of Eq. ([T]) for a single impurity interacting with electrons in an infinite 
volume. This is the traditional model for magnetic impurities in a metal. Note, however, that we ignore the presence 
of other impurities. A simplifying feature of this situation is that the mapping onto a single channel ID system is 
exact. As we will see in sub-section IIII Dt this is not generally the case when we put a 3D system into a finite box. 
Reflections off the boundary of the "box" greatly complicate the model. 

An immediate question that arises is: how dilute must the magnetic impurities be to justify a single impurity 
approximation? Given the screening cloud picture, one might come to the pessimistic conclusion that the average 
separation of impurities should be much larger than 

(Here n[ m ^ is the impurity density.) Such low densities are rarely, perhaps never, acheived in experiments on metals 
with dilute concentrations of magnetic impurities. Since typical ratios of £,k to lattice constant can be in the thousands 
this would require densities ?ii m p <C 10 -9 per unit cell. The reason 9 - that the condition is much weaker than this is 
basically that screening cloud wave-functions around different impurities become nearly orthogonal even when their 
centres are much closer together than their size. To see this consider two spherically symmetric wave- functions, centred 
around the points ±i?/2 but otherwise identical, only containing Fourier modes within a narrow band of wave- vectors 
of width around fcj?. The overlap is: 

0{R) = J d 3 rtp*(\f+R/2\)ijjQf-R/2\). (19) 

In terms of the Fourier transform of the screening cloud wave-function, this becomes: 

0<fl > - / !r ,£ *«*>i 2 = [ ikkIm)f w (20) 
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We expect the screening cloud wave-function to decay on the length scale Zk but to oscillate at the Fermi wave- vector. 
It should be built out of wave- vectors within Zk of k F . Thus we may assume: 

hK*OI 2 « f[(k - k F )u\ (2i) 

where the scaling function f(y) obeys the normalization conditon: 

J dyf(y) = 2tt 2 (22) 

in order that O(0) = 1. Thus, 

O(R) = (1/2tt 2 ) f dyf(y) sm[k F R + {R/£, K )y]/k F R. (23) 



For R <C £,Kj this reduces to: 

0{R) = sm(k F R)/k F R (24) 

independent of the details of the wave-function. Thus screening clouds centred on different impurities have negigible 
overlap provided that they are separated by a distance R^> kp 1 . A simple estimate of the condition for validity of 
the single impurity Kondo model is provided by the Nozieres exhaustion principle. There must be enough available 
conduction electron states, within of the Fermi wave- vector, to form linearly independent screening wave-functions 
around each impurity. This gives the conditions on the average impurity separation: 

» #V /3 m (25) 
n^ p 2 » ej 2 k F 1/2 (2D) (26) 
» Zk (ID). (27) 



A. Knight Shift 

The Knight shift is proportional to the r-dependent polarization of the electron spin density in linear response to 
an applied magnetic field. In NMR experiments, the total effective magnetic field felt by a nucleus is measured, from 
the nuclear resonance frequency. In addition to the applied magnetic field there is an additional contribution arising 
from the hypcrfinc interaction between the nuclear spin and the surrounding electron spins. Assuming for simplicity 
that this hyperfine interaction is very short range, the Knight-shift is just proportional to the local spin density at 
the location of the nucleus. For sufficiently weak applied fields, this is oc x( r )> the local susceptibility. In the absence, 
of the magnetic impurity, this is simply the r-independent Pauli susceptibility, xo- The Kondo interaction leads to 
an additional r-dependent term, which vanishes far from the impurity. Let 

5 d = J d 3 rS el (r), (28) 

be the total electron spin operator. The total conserved spin operator is: 

S tot = S imp + S e\- ( 29 ) 

The r-dependent magnetic susceptibility is: 

X (r) ee< S* d (r)S* tot > IT. (30) 

[We ignore the possible difference of (?-factors for impurity and conduction electrons for simplicity; this is discussed 
in Ref. [lOUllj .] x( r ) contains a constant term, present at zero Kondo coupling, which is simply the usual Pauli 
susceptibility, v/2. Since only the s-wave harmonic couples to the impurity spin, the other terms come entirely from 
the s-wave component of the electron field. It is convenient to take advantange of the boundary condition of Eq. (jlll) , 
to make an "unfolding transformation", regarding ipu(r) as the continuation of ipL( r ) to the negative r-axis: 



ipR(r) = ipL(-r). 



(31) 
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The s-wave part of the spin density operator may then be written: 

5fi wave (r) « ^[Jl(t) + Ji(-r) + e 2 * Fr 4W|^(-r) + +e- 2ikpr ^ L (-r)^ L (r)} (32) 

We expect this the continuum limit approximation to be valid at sufficiently long distances, r 3> kp 1 . Thus, we see 
that the local susceptibility takes the formH at r kp 1 : 



X {r) - u/2 -> - 1 [xun(r) + X 2fc F (r)2 cos 2fc F r] (33) 



Xun(r) = / rfr< [Vl(r,r)— Vi(r,r)+Vl(-r,r)— ^(-r,T)]5f ot (0) > 



where the uniform and 2fcp parts are slowly varying functions of r. These can be written in the ID field theory as: 

*r <r r»/»trr.^ 2 .. , , ... , 

X2kF (r) ee -v F J dr< 4(r,r) T ^(-r,T)5| ot (0) > . (34) 

Here S^ot is defined as in Eq. (1291) but we may restrict the electron part to its s-wave harmonic only, since the other 
harmonics give zero contribution to the Green's functions of Eq. (jM)) : 

4l=2^/ *l4(r)-V>L(r) + .... (35) 

It can be proved ) 10 ' 11 to all orders in perturbation theory, that 

Xun(r) = 0. (36) 

In the ID field theory, Xik F if) obeys a simple renormalization group equation with zero anomalous dimension which 
expresses how physical quantities change under a change of ultraviolet cut off and of bare coupling: 



Mi,(T,A,2>,rT/«f) = 0. (37) 



This simple form follows because only operators defined at r = have anomalous dimensions in a boundary field 
theory. Thus there is no anomalous dimension for ip L (r)ipL(~r). Furthermore, Stot has zero anomalous dimension 
because it is conserved, commuting with the full Hamiltonian. This RG equation implies that X2k F {T, Ao, D 1 rT/vp) 
which is, a priori a dimensionless function of 3 dimensionless variables, Ao, (the bare coupling), D/T and rT /vp, can 
be written as a function of the renormalized coupling at scale T, A(T) and rT /vf only. i.e. both the bare coupling, 
Ao and the ultraviolet cut-off, D can be eliminated in favour of a single variable, A(T). This is a basic consequence of 
renormalizability. An equivalent statement is that we may express Xik F as a function of the ratio tT/vf and A(r), 
the effective coupling at energy scale Vp/r. This follows from the RG equation, Eq. ^ which implies: 

F[X(r)} - F[X(T)} = \n(v F /rT) (38) 

where F(X) is the indefinite integral of 1 / /3(A): 

F{x) =rm dX '- (39) 

Eq. (|3"B"|) can be solved to express A(T) in terms of A(r) and a function of rT jvp- 

Results on X2fc F M were presented^ - - — to third order in pertubation theory in the Kondo coupling, Ao, at finite T. 
In the T — > limit these become: 

X2k F = ^I A o + A^ ln(Ar) + A^ ln 2 (Ar) + .5Ag ln(Ar) - Ag ln(X>/T) + constant • A§]. (40) 

Here A = AttD/vf- The presence of the A§ ln(D/T) term is very important. If this term were absent, we could safely 
take T — > without obtaining any infrared divergence in perturbation theory, at sufficiently small r. In that case, it 
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would be convenient to express X2k F in terms of A(r) only. At small A(r) perturbation theory would apparently be 
valid. However, the presence of this term implies that this doesn't work. We can write instead: 

X2k F « [A(r) + constant • A 3 (r)] [1 - A(T)] . (41) 

An important lesson from this expression is that the behaviour of the local susceptibity does not become perturbative, 
at low T, even at small r. This contradicts the naive expectation from the analogy with QCD, discussed in Sec. I. 
Nonetheless, it does suggest that there may be some crossover, or change in behaviour as we go from the region 
r <C £,k where A(r) < 1 to r > ^ where A(r) is large. A conjecture was made in Ref. [ll| for the behaviour at all 
r <C vp/T: 

X2kAr)^^f(r/U)x(T), (42) 

where x(T) is the total, r- independent, impurity susceptility; i.e. the linear response of the total magnetization to 
an applied field, with the bulk Pauli term subtracted off. Here / is some scaling function, depending on t/£k- This 
conjecture is consistent with both the perturbative results given above, Fermi liquid results at long distances and low 
T, r » T -C Tk , and also with results on the fc-channcl Kondo model in the large k limit. x{T) is well-known 
from a variety of methods, including BA. At T ^> Tk it has the pertubative behaviour: 



1 - A(T) 1 
47 1 4J" 



1 

1 



HT/Tk) 



(43) 



At low T < Tk, x(T) -> 1/(4T K ). For intermediate T it is quite well approximated by x{ T ) ~ -^/[{T + .65T K )]. 
We only know the form of the other scaling function, f(r/£x) asymptotically, at short and long distances. At short 
distances, r <C £k it is perturbative: 

f(r) - A W « _L_. (44) 

In the opposite limit, r 3> £if we may also determine^ f(r) from Fermi liquid theory. In the limit of low T <C Tk 
and large r>(jf first order perturbation theory in the Fermi liquid interaction of Eq. (1151) gives the proposed form, 
eq. (g2J) with: 

fir/be) -> 4r/«jc) + 1/3 + Ofor/r). (45) 

At intermediate length scales, of O(^k), the function f{r/^K) must somehow crossover between l/ln(r/^x) and 
4r/(7r£ft-). Determining the function in this region requires numerical methods. 

At T 3> Tr-, we found that the infrared divergences of perturbation theory are cut-off and weak coupling behaviour 
ensues, at any distance. There is no crossover at r ~ ^k in this case. Our result to third order in A is: 

X2k F {r) « (3tt 2 /4)A(T) 2 [1 - X(T)]e~^ rT ^ . (46) 

Our proposed scaling form for the local susceptibility when r <C vf/T: 

x (r)-u/2^^ff(r/U)x(T) (47) 

makes it clear that NMR measurements of the Knight shift would be a very difficult way of detecting the Kondo 
screening cloud. Typically an NMR signal is only picked up for nuclei within a few times kp 1 of an impurity spin. 
Over this range, the r-dependence of Eq. (|4"7|) is dominated by the first factor, cos 2k F r/r 3 . The slowly varying factor 
in the envelope of the oscillations, /(t/^k)) could not be measured unless signals could be picked up out to distances 
of 0(£if). At such large distances, the Knight shift is very small due to the factor of 1/r 3 . In fact, this form is 
qualitatively consistent with existing NMR experiments which looked for the Kondo screening cloudJ^ An important 
further complication in the experiments is that a given nucleus will feel a contribution to its Knight shift from many 
magnetic impurities at various distances. These contributions are expected to be simply additive, provided that the 
impurities are sufficiently diulte; see discussion at the beginning of Sec. |TTJ Nonetheless, it would be extremely 
difficult to disentangle these various contributions and require extremely accurate measurements of the Knight shift 
to determined the tiny contributions from distance impurities. 

An important lesson here is that the Kondo screening cloud is difficult to detect because it is so large. To see it, one 
must look for a crossover at large distances but at these distances the effects are weak precisely because the distances 
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are so large. The Kondo screening "cloud" is perhaps better described as a "faint fog" J£ It is worth noting that the 
situation is much better for a sample of reduced dimensionality. Our formula for the local susceptibility, Eq. (|33p . 
carries over directly to a 2D or ID sample, with the factor of 1/r 2 simply becoming proportional to l/r ^ 1 . In the 
ID case, x( r )/ cos(2kpr) crosses over from l/[rln(^/r)] for r < to a constant at r ^ £k, actually growing with 
distance. However, the diluteness condition is also modified in ID to 



"imp » & ( 48 ) 



as discussed at the beginning of Sec. HU 



B. Density Oscillations 



While the Kondo effect is associated with the spin degrees of freedom of the conduction electrons, it also has an 
interesting effect on the charge density in the vicinity of the impurityii^ This is perhaps a bit surprising, in light 
of spin-charge separation in ID. It turns out that density oscillations only occur when particle-hole (p-h) symmetry 
is broken. In particular, none occur in a half-filled tight-binding model, which has exact particle-hole symmetry. 
Importantly, the Kondo interaction itself preserves p-h symmetry. 

These density oscillations turn out to be quite a bit simpler to analyse than the Knight shift, discussed in 

sub-section III Al This is because, at long distances, r ^> kp 1 , they can be expressed as a half Fourier transform of a 
well-studied scaling function of w/Tk, the T-matrix, T(w). Due to the 5-function nature of the Kondo interaction, 
the exact retarded electron Green's function in the Kondo model can be written: 
i 

dte iut < t),^{f,0)} >=G{?,?,u) = G Q {r-r*,u) + G Q {r,u)T{w)G Q {-r*,u). (49) 



Here Go is the Green's function for the non-interacting case: 

f d 3 k 1 

G a (r,u>) = / 7 ^4e ife r — (50) 

J (27T)' j LU - 6 k + IT] 

The T-matrix is used to determined the resistivity and thus plays an important role in the theory of the Kondo model. 
Although it is not accessible by BA, NRG techniques have been developed to calculate it, in addition to perturbative 
and Fermi liquid results and other approximations. The density is obtained from the Green's function in the standard 
way: 

2 f° 

p(r) = / du)ImG(r, r, lo). (51) 

7T J-oo 

Note that G(f, r, ui) itself has trivial r-dependence, determined entirely by Go(r). Only the frequency-dependence 
reflects the Kondo physics. However, once we integrate over uj in Eq. (|51[) . this w-dependence introduces some 
interesting r-dependence. In the limit r ^> kp 1 , we may use the asymptotic expression for Go(r) 2 . In D-dimensions, 
for the free-electron dispersion relation this is: 

1 / — iki 



G (r, w)M 5- — — - exp(2ifc F r + 2iojrlv F ). (52) 

vp \ 2nr J 

This asymptotic form of Go can be calculated for other dispersion relations and has a similar form. Thus we see 
that, at r ^> kp 1 the density oscillations induced by a Kondo impurity are given by the half-Fourier transform of 
the T-matrix. Furthermore, for such large r, we expect this half Fourier transform to be dominated by frequencies 
< vp jr <C vpkp oc D. In this frequency range, the T-matrix is a universal scaling function of lo/Tk- 

It is important in calculating the density oscillations to take into account additional potential scattering, arising 
from the magnetic impurity, in addition to the Kondo interaction. This was ignored in sub-section III Al since it is not 
very important for the Knight shift. This corresponds to an additional term in the Kondo Hamiltonian, Eq. (|TJ) of 
the form: 

6H = Vj2i>U0)M0)- (53) 

a 

We refer to the s-wave phase shift at the Fermi surface due to the potential scattering as Sp. The T-matrix, at u> -C D 
can be expressed in terms of Sp, & universal function, txito/Tfc) and the density of states in D-dimensions, vp,: 

T(oj) = {e M - [t K (u/T K ) /{2-kv d ). (54) 
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Note that, in the particle-hole symmetric case, Sp = 0, T = tx / {2itvd)- By rescaling the frequency integration 
variable, and using v F /Tx — £,x > we see that the half- Fourier transform occuring in p(r) may be expressed in terms 
of a universal scaling function, F(r/£x)'- 

dwe 2iur / VF t K (u/T K ) = [v F /(2r)}[F(rT K /v F )-l}. (55) 

- oo 

Particle-hole symmetry of the model with no potential scattering (from which tx can be determined) implies that 
t* K (uj /Tx) — —tx{—w/Tx)- Furthermore, since tx is related to the retarded Green's function it is analytic in the 
upper half plane. These two conditions imply that F(r/£x) is real. On the other hand, note that it depends on both 
real and imaginary parts of tx- This makes it challenging to determine the density oscillations numerically, since 
while reasonably accurate results may exist for Im tx , this is not the case for the real part. From Eqs. (|5ip. (j5"4l and 
(|55[) we obtain the formula for the density oscillations at long distances from the impurity, r> kp 1 : 

p(r) - po ->■ ^-[cos(2k F r - ttD/2 + 26 P )F(r/£ K ) - cos(2fc F r - txD/2)] (56) 

where the constant is C3 = l/(47r 2 ), C2 = l/(2ir 2 ) and C\ = l/(27r). The formula applies in D=l,2 or 3, with a 
spherically symmetric dispersion relation. It also applies to the ID tight-binding model in which case r is restricted 
to integer values. Note that in this case, there is exact p-h symmetry at half-filling with no potential scattering, 
corresponding to kp = tt/2 and bp = 0. In this case, the right hand side of Eq. (|56[) vanishes identically as expected, 
p-h symmetry breaking is necessary to get non-zero density oscillations and the potential scattering phase shift, Sp 
plays an important role in them. 

The universal Kondo T-matrix, tx(u /Tx), can be calculated at higher frequencies, uj ^> Tx using perturbation 
theory in the Kondo interaction. It can be calculated at low frequencies using Fermi liquid theory. We find that these 
expansions can be simply inserted into Eq. (|55[) to obtain expansions for F(r/£x) and hence the density oscillations 
at short distances, r <C £x and long distances, r>(j;. The weak coupling result is 

t K (u) = -(3iir 2 /8)[\ 2 Q + 2X 3 Q H V /u) +•••]■ ( 57 ) 

This may be written: 

txiu/Tx) -> -(3ttV8)A 2 (M) a -to?i/[Z\£{\u\/T K )]. (58) 

The T-matrix vanishes slowly at high frequencies. Substituting this weak coupling expansion into the integral defining 
F, Eq. (1551) and integrating term by term gives: 

F = 1 - (3tt 2 /8)[A^ + 2Xl ln(Ar) + . . .] (59) 

where A is of order v F /D. We recognize the weak couplig expansion of the effective coupling at scale r: 

F(r/£ K ) -> 1 - (3vr 2 /8)A 2 (r) « 1 - 3vr 2 /[8 ln 2 (^/r)]. (60) 



Note that in the weak coupling limit F — > 1 and then Eq. (|56[) reduces to the standard formula for the Friedel 
oscillations induced by a non-magnetic impurity which produces a phase shift at the Fermi surface of Sp (in the 
s-wave channel only). 

Our results indicate that weak coupling behaviour occurs for the density oscillations at short distances r <C £,x even 
at T — 0. This is quite different than what we found for the Knight shift in Sec. Ill Al In that case we found, in Eq. 
(|40p. that the infrared divergences of perturbation theory were not completely cut off at short distances. However our 
conjecture, Eq. (|42j) . implied that the non-perturbative aspect simply led, at low T -C v F /r, to a factor of x(T) in 
X(r) which was otherwise linear in A(r), Eq. (j44|) . The density oscillations therefore show simpler behaviour at short 
distances. 

At low frequencies, uj <C Tx, the Kondo T-matrix can be calculated from Fermi liquid theory: 

t K -> -*[2 + iTTLu/(2T K ) - 37r 2 w 2 /(I6T|.) + ...]. (61) 

Substituting into Eq. (1551 . gives: 

F(r/£ K ) -► -1 + Trfr /(4r) - 37r 2 ^/(32r 2 ) + . . . (62) 

At large distances, r ^> £xi F — > — 1. From Eq. (1551) . we see that this is equivalent to shifting Sp — >• S P + tt/2 while 
setting F = 1, the non-interacting value. In this limit the Kondo scattering simply changes the phase shift by 7r/2, 
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the familiar result. Thus, while we obtain the standard Friedel oscillation formula at both short and long distances, 
there is a non-trivial and universal crossover, determined by F{r/^x), at intermediate distances. 

The density oscillations were calculated using an improved NRG method in Ref. [7] . In this approach Wannier states 
are introduced both at the impurity location and at the point of interest, r, fixing the problem with spatial resolution 
in the usual NRG approach. Results are shown in Fig. (Q}. The oscillations can be parametrized, at r ^> kp 1 , by 
the form of Eq. (f56|) . The function F(t/^k) fits the asymptotic predictions at large and small arguments, crossing 
over between 1 and —1 as expected. In fact, it is fairly well fit, throughout the crossover by the simple "one spinon 
approximation" 16 ' 17 



t K (u/T K ) « -2i/[l - m/(2T K )} 

giving an approximate formula for F: 

F(u) sa 1 + (8u/7r)e 4 " /,r Ei(-2u/7r), 
where Ei(x) is the exponential-integral function. 



(63) 



(64) 
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FIG. 1: NRG results on charge oscillations around a Kondo impurity coupled to ID conduction electrons with particle-hole 
symmetry from Ref. 7|. Note that the oscillations vanish at kFr/n £ N. As shown in the inset, the properly rescaled envelope 
function of the oscillations (extracted as p — po at the local maxima) for different Kondo couplings nicely collapse into one 
universal curve except for the points where r ~ kp 1 . In the inset we show the analytical results for the asymptotics as well: 
Note the good agreement between the analytical results and the numerics. 

Experimental measurement of the Kondo screening cloud via density oscillations would again be extremely chal- 
lenging. The l/r D factor in Eq. (f5fj|) implies that the oscillations would be extremely small at the length scale ^k- 
(Note that the number of oscillations that would need to be measured is of order ^k^f ^> 1-) As for the Knight shift, 
the situation is improved in lower dimensions. The Kondo effect is apparently observed for magnetic impurities on 
metallic surfaces using Scanning Tunelling Microscopy (STM). In cases, where the Kondo interaction is predominantly 
with surface conduction electron states, Eq. (|56l) with D = 2 should apply. The STM tunelling rate for tip energy E 
is usually assumed to be proportional Eq. (|51l) but with the lower limit of integration replaced by E. This quantity 
could also be calculated from knowledge of , but only when E ^> Tk do we recover a scaling function of rj £,k only. 



C. Impurity spin correlation function 



Another quantity which gives a very direct picture of the Kondo screening cloud is the equal time ground state 
correlator of the impurity spin with the spin density at a distance r: 



K(f,T) =< 0\S* el (f,0)S? mp (0)|0)> 



(65) 
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Although this quantity was also discussed at finite T in Ref . [lolflTj , we restrict our discussion here to the T = case. 
Unlike for the local susceptibility, no divergences are encountered in perturbation theory by taking T — > 0, at least to 
the (third) order studied. K(r) can easily be seen to obey an exact sum rule. If we consider a finite system with an 
even number of electrons in total (including the impurity spin as an electron) then the ground state is a spin singlet: 

Stot(0)*|0>=0. (66) 

Thus we have: 

=< 0|Sf mp (0)S tot (0)*|0 >= 1/4 + J d\K{r). (67) 
A heuristic picture of K{r) is obtained by writing: 

S^if) « ^p(r) (68) 

Here p(r) is the probability of the screening electron being at the location f and s^j is a S=l/2 operator, representing 
the spin of the screening electon, with < S? m >= —1/4. p(r) = \<fr(r)\ 2 where 4>(r) is the screening cloud wave- 
function. Then we obtain K(r) s» —p(r)/A and the sum rule of Eq. (|67p is simply the normalization condition on the 
screening electron's wave- function. As for the local susceptibility, discussed in sub-section III A[ at r ^> kp 1 we may 
decompose K into uniform and 2kp parts: 

K ( r ) o 2 1 2 [KMr) + 2K 2kF (r)] (69) 

where -ftTun and K2k F can be calculated in the ID field theory. Perturbation theory in the Kondo coupling, valid at 
distances r <C £k, gives : 10 ! 11 

v 7r^A(r) 2 (l + A /2) 

Kun 2r 

_ ^ F A(r)(l + A /2) 
K2k F ^ (70) 

Several differences are to be noted with the local susceptibility. First of all, we have obtained a finite result at T = 0, 
as noted above. Secondly the uniform part is now non-zero. Thirdly, the result cannot be expressed in terms of 
the renormalized coupling at scale r only, but also involves a factor containing the bare coupling, Ao- This can be 
understoo d 10 ! 11 from the fact that the impurity spin operator has a non-zero anomalous dimension: 

7imp «A 2 /2 + ... (71) 

However, this is not very important since this factor goes to 1 in the scaling limit of weak bare coupling. On the other 
hand, at r > we may calculat o 10 ' n i 18 K{r) using Fermi liquid theory giving: 

K 2kF -> -(l/2)tfun const x (72) 

Our heuristic interpretation then suggests a slow power-law decay for the the screening cloud electron probability 
density: 

/ \ sin 2 kprfx , * n 
p(r) « const x r » U . 73 

On the other hand, at shorter distances, r <C K(r) is not negative definite and the heuristic picture of K{r) is 
not valid. 



III. SPIN IN A MESOSCOPIC DEVICE 



Given the extreme difficulty of observing the Kondo screening cloud from length-dependent measurements in a 
macroscopic sample, discussed in Sec. [Hi together with the fact that mesoscopic devices routinely contain components 
with dimensions of order 1 micron, it is natural to consider whether the best way of finally observing the Kondo 
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screening cloud experimentally might be via size-dependent effects in such a device. We are interested in devices 
containing a quantum dot in the Coulomb blockade regime with a gate voltage tuned so that the number of electrons 
in the dot is odd, with a S=l/2 ground state. For simplicity, we will also assume that the tunelling from dot to leads 
is sufficiently weak that only virtual charge fluctuations of the dot need be considered so that the Kondo model is 
appropriate, rather than the Anderson model. Alternatively, we might consider one or more magnetic impurity atoms 
in a mesoscopic sample. 

An important point is that, associated with a finite sample size, we have a finite spacing between energy levels. If 
the sample is 1 dimensional, then the level spacing will generally be A vf/L where L is the sample size. Thus the 
condition vf/Tk w L is equivalent to A w Tk- Thus, it sometimes argued that observing size dependence in this 
situation does not really show the existence of a Kondo screening cloud but "merely" that Tjc is the characteristic 
energy scale for the Kondo model. However, since as discussed in the previous sections, the Kondo screening cloud 
concept is really just the inevitable consequence of being able to convert an energy scale to a length scale using a 
factor of velocity, size effects can provide realizations of it. We may equally well say that the behaviour changes when 
the finite size gap becomes larger than Tk or else when the Kondo screening cloud no longer fits inside the sample. If 
the sample should instead be regarded as 2 or 3 dimensional the finite size level spacing may be <C vf/L depending 
on which energy levels are important and other details, discussed in Sec. IIIID1 

The Kondo model is used both to describe magnetic impurities in metals and also gated semi-conductor het- 
erostructure quantum dots. The heterostructure (such as GaAs/AlGaAs) provides a fairly clean 2D electron gas 
(2DEG) buried some distance (often around 100 Angstroms) below the surface of the semi-conductor wafer. Gate 
voltages are applied to the surface to define point contacts and quantum dots. A quantum dot refers to an island of 
electrons, with the electron number (typically around 100 or less) controlled in unit steps by a gate voltage. When 
this number is odd, the quantum dot usually has an S=l/2 ground state. The quantum dot can be connected by 
narrow point contacts to the left and right sides of the 2DEG. If the point contacts are close to being pinched off they 
only permit one channel of electrons to pass through, giving a 2e 2 /h conductance. In this case, a simplified model of 
the system involves effectively ID leads. The quantum dot has a charging energy, U, as well as a gate voltage, such 
that the energy as a function of electron number, N, is U(N — No) 2 for a value of the parameter Nq controlled by 
the gate. It is often convenient, theoretically, to use a ID tight-binding model; universality of Kondo physics implies 
that such details are not important. The corresponding Hamiltonian is: 

H = H Q + H d (74) 

where 

H = - 

and 

H d = eo"o + U n Qt n Qi . (76) 

Here riQ a is the occupation number at site for spin a, with no = not + n oi- t' represents the tunelling amplitudes, 
through the left and right point contacts, to the leads. (We take these to be equal for simplicity but the generalization 
is straightforward.) 

If the tunnelling amplitudes are sufficiently weak compared to U, then the quantum dot makes only virtual charge 
fluctuations and the system is in the Coulomb blockade regime. The conductance through the dot then usually 
begins to decrease as the temperature is lowered, if No is close to an integer. In the case where No is near an odd 
integer, so that the quantum dot behaves as an S=l/2 impurity, we may make a Shricffcr- Wolff transformation to the 
corresponding Kondo model with Hamitonian: 

H = H a + H K (77) 

where now: 

(78) 



j=— oo 



l,a - t C ja C j+l,a> ~ *' ' c -la c 0a ~ t ' 4) a C la + h.C. 
J'=l 



(75) 



Ho 



-2 

< E 

j=—oo 



H K = J(C-l + C^-id + C_i) • S 



(79) 
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where S is the spin operator for the site 0: 



t CT 



and the Kondo coupling is: 



J = 2r 



1 



<7 + e 



(80) 



(81) 



In general, a potential scattering term is also generated, of the same order of magnitude as J. To simplify the 
discussion, we focus here on the particle-hole symmetric case, eo = —U/2, where the potential scattering vanishes. We 
also assume that the system is at half-filling so that the particle-hole symmetry is exact. It is then easy to calculate 
the zero temperature conductance for this Kondo model in either limit, J <C t or J 3> t. The second limit is not 
physical for the underlying Anderson model of Eq. (|74[) but it is convenient to consider nonetheless as a low energy 
fixed point Hamiltonian. When J = the the Kondo Hamiltonian contains no terms linking left and right sides of the 
dot, so the conductance vanishes. When J>t, one electron gets trapped in the symmetric orbital, with annihilation 
operator: 



c s = (ci + c_i)/V2 



(82) 



and forms a singlet with the impurity spin. At first sight, one might think that this would block the transport through 
the impurity. However, this is not so due to transport through the antisymmetric orbital with annihilation operator: 



c a = (ci - c_ 



(83) 



Electrons entering from sites —2 or 2 can hop into this antisymmetric orbital without breaking the Kondo singlet. If 
instead they hop into the symmetric orbital there is a large energy cost of 0{J). At large J/t we may obtain a low 
energy effective Hamiltonian from Hq by projecting c±\ onto c Q : 



Pc ±1 P = ±c a /V2. 



This gives: 



Hi 



low 



-t 



-3 

E 

j=—oo 



C ia C j+l,<* 



j=4 



71' 



.C. 



The tranmission probability for this non-interacting electron problem is easily calculated: 20 

T(k) = sin 2 k. 



(84) 



(85) 



(86) 



In particular, at half-filling, kp — 7r/2, the transmission probability at the Fermi energy is one (corresponding to a 
resonance). The conductance for this non- interacting model is given by the Landauer formula: 



G 



2<r 



T(e F ) 



2e* 
h ' 



(87) 



Calculations show that, for a small bare Kondo coupling, the conductance grows with decreasing T, saturating at 
the perfect tranmission value, 2e 2 /h, for T -C Tk- These results can be obtained by expressing the conductance as a 
frequency integral of the imaginary part of the T-matrix, introduced in sub-section III Bt at finite temperature, times 
the derivative of the Fermi function. Analytic formulas can be obtained for G(T) in both the high T and low T limits 
using perturbation theory and Fermi liquid theory respectively. Calculations at intermediate T are typically based on 
less controlled approximations. This increase of conductance with decreasing T, due to the Kondo effect, only sets in 
at low temperatures, following an initial decrease due to the onset of the Coulomb blockade. 

Note that this behaviour is the inverse of what happens for a magnetic impurity in a bulk metal where it is the 
resistivity which grows with lowering T, not the conductance. The impurity behavior is more closely related to that 
for a side-coupled quantum dot where the conductance is 2e 2 //i for zero Kondo coupling (high T) and vanishes at low 
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A. Persistent current in a ring containing a Kondo impurity 



A calculationally simple situation in which to observe finite size effects is to close the leads embedding the quantum 
dot into a ring. The Hamiltonian is that of Eqs. (|77H79I) but the leads are now of finite length with periodic boundary 
conditions: 

L-2 

H = -t^2(c] a c j+lta + h.c.) (88) 

H K = J(c L _i + ci)t|( Ci _i + ci)-^ (89) 

While the conductance is now not measurable, one can instead calculate the persistent current in response to an 
enclosed magnetic flux, $ = (Hc/e)a. The current, at T = 0, is given by the derivative of the ground state energy 
with respect to the flux: 

j = -(e/h)dE /da. (90) 

A perturbative calculation gives: 19 ' 20 

j e (a) = ^^{[sin5[A + A 2 ln(Lc)] + (l/4 + ln2)A 2 sin25} + 0(A 3 ) 

3o(a) = ^^sin2a[A 2 + 2A 3 ln(Lc')]+0(A 4 ), (91) 
16L 

for A even and odd respectively, where c and d are constants of 0(1) which we have not determined and: 

a = a (A/2 even) 

= a + n (A/2 odd). (92) 



O: 



A is the number of electrons, including the electron on the quantum dot, j — 0. At half-filling, A is just the total 
number of lattice sites, including the origin; N = L. The fact that j is 0(A) for A even but not for A odd is easily 
understood. The unperturbed ground state consists of a partially filled Fermi sea and a decoupled impurity spin. For 
A even, there are an odd number, A — 1, of electrons in the Fermi sea. The unpaired electron at the Fermi surface 
forms a spin singlet with the impurity in first order degenerate perturbation theory. On the other hand, for A odd, 
there are no unpaired electrons in the non-interacting Fermi sea so it is necessary to go to second order in A. A very 
important property of these results is that, to the order worked, the persistent current only depends on the effective 
Kondo coupling at the length scale L: A(L) = A + A 2 InL. That is, logarithmic divergences are encountered in next 
to leading order, as is standard for many calculations in the Kondo model, but they only involve InL, not InT. The 
current is finite at T — 0, for finite L. This suggests that the finite size of the ring is acting as an infrared cut-off 
on the growth of the Kondo coupling. Provided that A InL <C 1, the higher order corrections are relatively small 
and perturbation theory appears trustworthy. This is equivalent to the condition 3> L. We may say that, in this 
case, the screening cloud doesn't "fit" inside the ring so the Kondo effect (growth of the effective coupling to large 
values) doesn't occur. We expect that higher order perturbation theory would preserve this property, giving Vpe/L 
times functions of the renormalized coupling at scale L (and the flux) only. This follows because the current obeys a 
renormalization group equation with zero anomalous dimension. This in turn follows from the fact that the current 
is conserved d < j(x) > jdx — so it can be calculated at a point far from the impurity spin. This scaling form 
implies, equivalently, that we can write the current as: 

Je/o = ^fe/o{L/£,K,&) (93) 

for even and odd A respectively, where f e and f a are universal scaling functions, depending only on the ratio L/^k 
and not separately on the bare Kondo coupling Ao and cut-off T>. In the limit L ^> the persistent current can be 
calculated using the non-interacting effective low energy Hamiltonian of Eq. i.e. Fermi liquid theory. This gives: 

j e (a) = -^£[a-7r] (A even), (94) 
j a (a) = + [ 7r - 7 r])] (A odd). (95) 

7rL 



15 



Here 

[a] = a (mod 2tt), \[a}\ < n. (96) 

These small L and large L limits are plotted in Fig. ([2]). Note that the current has the same sign for all a in both 
limits. Furthermore it has the same period: 2tt for N even and 7r for N odd. For intermediate lengths, L of order 
it is necessary to do a numerical calculation. A combination of exact diagonalization and DMRG [21| supports the 
scaling behavior of Eq. (|93|) giving scaling functions that interpolate smoothly, as a function of £k/L. See. Fig. ((3}. 




-2 



FIG. 2: Persistent current vs. flux for an even or odd number of electrons in the weak coupling limit, for £k/L 50, from Eq. 
(|91fl (solid line) and in the strong coupling limit, £k/L « 1, from Eq. (|95[) (dashed line). j is multiplied x5 at £ k/L — 50 
for visibility. The solid lines are obtained from Eq. (|9ip using the effective coupling X(L) « 1/ \xi{^k/L). (From [IJa].) 




FIG. 3: The current at 1/2-filling for iV = 4p — 1 (A) and N = Ap (B) for the embedded quantum dot. Results^ are shown 
for a number of system sizes with Jk = 1 as a function of a/n. The results for the small system sizes are obtained using exact 
diagonalization methods while the larger system sizes (circles) have been obtained using DMRG techniques. From Ref. [2lj . 
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B. Charge steps in a finite length quantum wire terminated by a quantum dot 

Another mesoscopic system which is quite readily analysed theoretically involves a linear quantum wire of finite 
length with a quatum dot, in the Kondo regime, at one end. [See Fig. Q.] We assume that the system is very weakly 
tunnel-coupled to an electron reservoir and that a uniform gate voltage, Vq, is applied to the wire, corresponding to a 
chemical potential /i = —eVa, in addition to the gates controlling the occupancy of the quantum dot and the tunelling 
between dot and wire. As fj, is varied, the number of electrons in the system, in its ground state will vary in single 
steps at particular values of fi. This defines a "charge staircase". See Fig. ([5]). The width of the steps, with an even 
or odd number of electrons in the system, is a sensitive indicator of the strength of the effective Kondo coupling. The 
steps, at discrete values of n, could be detected by some Coulomb blockade technique. Such measurements are likely 
to be far easier than measuring persistent currents in closed rings. For simplicity, we treat the wire as having only a 
single channel - i.e. as being an ideal ID system. A great advantage of this system is that the charge staircase can be 
determined very accurately from the Bethe ansatz solution of the Kondo model. Remarkably, this remains true, 
to some extent, even when we include short range Coulomb interactions throughout the wire. In this case we take 
advantage of the more recent Bethe ansatz solution of an S=l/2 Heisenberg antiferromagnet with a weakly coupled 
spin at one end (22[. We first discuss the case with no Coulomb interactions (apart from those in the quantum dot, 
leading to the effective Kondo model) and later include bulk Coulomb interactions in the wire. 

k l >| 



v du 



FIG. 4: Possible experimental setup. Vdw 
varies the chemical potential in the wire. 
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FIG. 5: Charge quantization steps for the wire coupled to a small quantum dot. The arrows indicate the direction in which 
the single steps move as A (L) grows. 
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1. Interactions in quantum dot only 

Again it is convenient to consider a tight-binding model. In the Kondo limit this is 

H = H a + H K (97) 

with 

L-2 

H = -tJ2 [(c]« c i+i.« + h - c -) - ^ c ] c i\ ( 98 ) 

and 

H K = Jc\^ Cl -S. (99) 

Note that this time we have "open" boundary conditions. The chain terminates at site L — 1. We are interested in 
how the total electron number in the system, N, changes at T — in the grand canonical ensemble, as we vary /j,. 
We define N to include the single electron on the quantum dot. Considering a small range of fi, the steps with even 
N all have the same width, as do the steps with odd TV. We are interested in the ratio R: 

Rm'g (100) 

where i5/i e /° is the width of the interval of chemical potential over which N has a fixed even (odd) value. 

The limits of weak and strong Kondo coupling are readily analysed^ See Fig. ([5]). When J — 0, finite width steps 
only occur for N odd since the energy levels in the quantum wire are doubly occupied and a single electron resides 
on the quantum dot. In the opposite limit of strong Kondo coupling, finite width steps only occur for N even. One 
electron is removed from the Fermi sea to screen the spin and the remaining electrons doubly occup free fermion 
energy levels, appropriately tt/2 phase shifted. Thus R — oo at zero coupling and R = at strong coupling. Again it 
is possible to calculate the corrections to these limits at leading order in renormalization group improved perturbation 
theory at weak coupling and Fermi liquid theory at strong coupling. At weak coupling we find: 

|»(3/2)[Ao + Agln(L/c)] + ... (101) 

As for the persistent current, we see that the length of the quantum wire cuts off the growth of the Kondo coupling, 
at T = and we may write: 

1 -> (3/2)A(L) w - 3 (L«Ck)- (102) 

R 2m(£ K /L) 

In the other limit, first order perturbation theory in the Fermi liquid interaction of Eq. (1151) gives: 

R ^^T' (103) 

Again, we see that R appears to be a scaling function of £k /L only. 

For this case, we can check this scaling hypothesis, and determined the function R(L/£k) to high precision for all 
L/£k by using the Bethe ansatz solution4£ The result is shown in Fig. ©. Excellent agreement is obtained with 
the asymptotic weak (L <C £,k) and strong <C L) coupling limits. Weak particle-hole symmetry breaking has a 
negligible effect on these results, since it simply shifts the entire charge staircase by a constant^ 



2. Including Coulomb interactions in the quantum wire 



So far, we have ignored the effect of bulk Coulomb interactions, only considering the Coulomb interaction on the 
quantum dot, in the Anderson impurity model, which leads to the Kondo interaction. This neglect may be justified 
in many cases, based on Landau Fermi liquid theory and the fact that the Kondo interaction just involves low energy 
electron states. However, it is certainly not justified when the Kondo screening takes place in a ID quantum wire. 
In this case, Fermi liquid theory breaks down, being replaced by Tomonaga-Luttinger Liquid (TLL) theory. In ID, 
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FIG. 6: Universal ratio of even and odd charge steps, R — S/j,°/Sjj, c , as a scaling function of £k/L from [2^] - Bethe ansatz 
results-- obtained for various systems sizes (No = 51, 101, 201, 501, 1001, 2001) and 13 different values of the Kondo exchange 
J, indicated by different symbols. For each value of J, the system lengths L have been rescaled L — > L/£k(J) in order to 
obtain the best collapse of the data using the strong coupling curve Eq. (|103p (dashed red line) as a support for the rest of 
the collapse. The weak coupling regime for L <C £k, enlarged in the inset (a), is described by the weak coupling expansion 
Eq. (|101|) with constant ~ 0.33 (continuous blue curve). Inset (b): The Kondo length scale, extracted from the universal data 
collapse of the main panel (black squares), is described by theexpected exponential dependence on Kondo coupling (dashed 
green line). 



Coulomb interactions lead to major changes in the interaction with an impurity. It turns out however, that these 
changes are minimized when the magnetic impurity is at the end of a quantum wire, so that the conclusions of the 
previous sub-sub-section will not be modified too drastically. 

We modify the model of Eq. and (IM1) by adding screened bulk Coulomb interactions. For example, we could 
consider the Hubbard model with only on-site repulsion: 

H c = U^2n jt n ji (104) 

3 

where fij a is the number operator for electrons of spin a at site j. We expect our results to apply to more realistic 
models, not based on tight-binding approximations and with longer range screened Coulomb interactions. We begin 
by writing the ID electron annihilation operators in terms of left and right movers, similar to what we did for the 
s-wave projected 3D theory in Eq. (|T!Zj) . 

Cj « e ik ^M3) + e~ lkFJ Mj)- (105) 

Here iPr/l vary slowly on the lattice scale, containing Fourier modes of the Cj's in narrow bands around zt&F, 
appropriate for the low energy effective theory. The Coulomb interactions reduces, in this low energy approximation, 
to four different terms which can be conveniently written in terms of charge and spin densities (or currents) for left 
and right movers: 

jR/L = ^2^ R/La ^R/L a , JR/L=^2^ R/La ^-^R/L0- (106) 

a a/3 

Three of the Coulomb interaction terms are proportional to (J\ + J%), JlJr and ( j£ + J%)- To make further progress, 
we bosonize the fermions, and introduce separate bosons for the charge and spin degrees of freedom. Upon doing 
this, the first two interaction terms only involve the charge boson. They change the velocity, v c for charge excitations 
and also rescale the boson field by the Luttinger parameter for the charge sector, K c . The third term changes the 
velocity, v s for spin excitations. Importantly, all of these terms leave a non-interacting model of two massless bosons, 
for charge and spin. The fourth Coulomb interaction (exchange) term is written explicitly as: 



H ex = ~2ng v s / dxf L ■ Jr, 



(107) 
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where the coupling constant, go, is proportional to the strength of the Coulomb interactions. This can be written 
entirely in terms of the spin boson operators but corresponds to a non-trivial interaction. The exchange interaction, 
go is marginally irrelevant, its effects getting weaker at lower energy, so that asymptotically, free boson behaviour is 
obtained. Note that this is the opposite to the behaviour of the Kondo coupling which gets stronger at lower energies. 
We assume the system is not at half-filling so that we can ignore Umklapp interactions. 

An important consideration for the charge steps in the wire-dot system we are considering is the boundary conditions 
on the low energy fields, 4>r/l at j = 1, the end of the chain next to the quantum dot. In the limit, J = these are 
simply: 

V>l(0)+Va(0) = 0. (108) 

These can be seen, for example, by adding a "phantom site" at j = and requiring that c = 0. Since the low energy 
fields ipR/L var Y slowly over one lattice spacing, the Kondo interaction can be written in terms of Jr(0) only: 



ci « *p R (0y kF +M0)e- lkF «2isinfcWfl(0) (109) 
ci w 4sin 2 fc F Jij(0). (110) 



Since Jr can be expressed in terms of the spin boson only, the Kondo interaction only involves the spin degrees of 
freedom. This is not the case if the impurity interacts with the quantum wire far from its end. In that case, we must 
use the bosonization expression: 

«) « e 2ikFi exp{i^K~ c <b c ]n (111) 

where n can be expressed in terms of the spin boson. Thus, in this situation the Kondo interaction mixes the charge 
and spin sectors leading to quite different physics. Here, we have deliberately considered the case of the end-coupled 
spin in order to preserve as much as the physics of the usual Kondo model as possible. 

The charge staircase is determined by the ground state energy, as a function of number of particles. In the case 
where we set to zero both the Kondo interaction, Ao and the bulk marginal coupling, go, this is: 

E Q -^N=^ Z (N-N r + ^ (112) 

where s is the spin quantum number of the ground state: for N odd and 1/2 for N even. (Recall that we include 
the electron on the decoupled quantum dot in our definition of N. Here s is the spin of the electrons in the wire, not 
including the quantum dot.) Including both exchange and Kondo interactions, it is convenient to write: 

E -fxN= -^(N - N ) 2 + ^ /s (A ,.go,£) (113) 

where s — or 1/2 for N odd or even and f s are some functions to be determined. Then defining: 

/ = A/2 - fo (114) 



it is easily seen that the ratio of step widths is: 



where 



*-l±3 (IK) 



u = v s K c /v c . (116) 

Coulomb interactions made v c < v s rts v F and also reduce K c from its non- interacting value of 1. Both these effects 
reduce u from its non- interacting value of 1 thus bringing R closer to 1. For cleaved edge overgrowth quantum wires 
MRi.5 and for carbon nanotubes u w .1. 
In lowest order perturbation theory, 

/ « 1 - lg(L) - 3A(L), (£«6f) (117) 
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Here we have replaced both the exchange and Kondo couplings by their renormalized values at scale L. The exchange 
interaction lowers the ground state energy when the wire (excluding the quantum dot) has an s = 1/2 ground state, 
corresponding to N even. The effective exchange coupling is: 

s(L)«l/la(L/Li) (118) 



where L\ is a short length scale, of order k F x or smaller. The exchange interactions actually modify the /3-function 
for the Kondo coupling: 

X 2 +g\ + ... (119) 



d\nL 

(On the other hand, the Kondo interaction doesn't affect the rcnormalization of g(L) since the Kondo interaction 
only occurs at the end of the wire and the exchange interaction occurs throughout the wire.) Note that the exchange 
interaction causes the Kondo coupling to grow faster as the length scale is increased and consequently reduces £k for 
a given bare Kondo coupling: 



U oc exp[v/a/A + b] (120) 

for some constants a and b depending on g . For long lengths, L>(x, we may again use Fermi liquid theory. Note 
that the effective exchange coupling, g(L) continues to get smaller as we increase L so we may continue to treat it 
perturbatively. When however, one electron is removed from the Fermi sea to screen the spin of the quantum 

dot so that the remaining nearly free fermions system has one fewer electron. This switches the s = and s = 1/2 
ground states and changes the sign of the term in / linear in g(L): 

f^-l + L(L) + ^f, (121) 



In general we may write: 



2 JV ' 2L 



f = f(U/L >g (L)). (122) 



The dependence on g(L) destroys the universal scaling with L/£k but in practice this effect can be small since g{L) 
tends to be small and slowly varying. 

Since the Kondo effect for a Luttinger liquid, with the impurity spin end-coupled, only involves the spin degrees of 
freedom, it is essentially unchanged if the charge degrees of freedom are gapped, as would occur in the Hubbard model 
at half-filling. In fact, we may take the large U limit at half-filling to obtain the S=l/2 Heisenberg antiferromagnetic 
chain with one weakly coupled impurity spin at the end. The only important effect of this large U limit is that the 
bare exchange coupling, go becomes 0(1). This impurity model has been solved via Bethe ansatz methods^ and 
this allows an exact calculation^ of the function f(L/£jc,g(L)) for a particular (large) go. The results are shown, 
for L = 1000 and a range of Kondo couplings, in Fig. ([7]). The expected asymptotic limits, ±(1 — 3g(1000)/2) are 
obtained for L/£k "C 1 and 3> 1. Note that, even for this large bare exchange coupling, the deviations from the model 
with no bulk Coulomb interactions is quite small. For conducting quantum wires, we expect go to be considerably 
smaller and the deviations to be even less. Thus the main effect of the Coulomb interactions is to reduced R due 
to the parameter u in Eq. ([1151) . If R is the step width ratio in the non- interacting model, we may write it in the 
interacting case as approximately: 

R{U,L) ~ (123) 

This varies from R(£k/ ' L) = R{^k/L) at u = 1 (no Coulomb interactions) to R(£,k/L) = 1 at u = (strong Coulomb 
interactions). Varying £,k/L at fixed u, R goes from (l + u)/(l — u) at L <C £k to (1 — u)/(l + u) at 3> L. Coulomb 
interactions weaken the even-odd alternation in the step widths but only trivially modify the dependence on ^k/L. 



C. Transport through a finite wire containing a quantum dot 

The previous examples have involved closed systems, making the measurement of the Kondo effect challenging. Most 
experiments observing the Kondo effect in quantum dots involve transport measurements; the systems are "open" , 
connected to external leads. A natural way of observing finite size effects in a transport measurement is to embed a 
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FIG. 7: Effect of the marginally irrelevant bulk operator on the scaling function f(£i</L,g(L)) calculated from the Bethe 
ansatz from [23]. The symbols correspond to the noninteracting case, g{L) = 0, as in Fig. [6] The black circles are obtained^ 
with the Bethe ansatz solution of the spin chain Kondo model for L — 1000 (<?(1000) ~ 0.115). 



quantum dot, in the Kondo regime, in a finite length quantum wire which is weakly connected to external leads by 
tunnel junctions ) 24 ' 25 as illustrated in Fig. (JH). Assuming for simplicity, perfect parity symmetry, the Hamiltonian is: 



H — Hi + H w + H LW + H K 
where L labels the leads, W labels the quantum wire and K the Kondo interaction. Here: 



H, 



-t 



-L-2 oo 

E + E 

j=-oo j=L + l 



-2 L-l 

11 = 1 E + E 

j=-L j=l 



(c]c j+ i + h.c.) 
(c]cj + i + h.c.) + e\ 



-1 L 

E+E 

j=-L j=l 



Hlw 



(124) 

(125) 

(126) 

(127) 
(128) 



Here spin indices are suppressed and we have ignored a potential scattering term, as well as Coulomb interactions in 
the wire, for simplicity. Note that we have included a chemical potential term, e\y for the wire. When tzw = the 
wire + dot system is closed and we expect that the length, L, of the wire will act as an infrared cut-off on the growth 
of the Kondo coupling, leading to scaling in £,k/L as discussed in the previous two sub-sections. What happens when 
the wire is connected to the external leads? When <C L, the Kondo screening cloud fits inside the leads and, at low 
T <C Tk, we therefore expect perfect transmission through the quantum dot. The conductance is then simply that of 
a finite length non-interacting wire, W, connected by tunelling amplitude twh to the external leads. The conductance 
is then easily calculated by the Landauer formula which expresses it in terms of the transmission probability for an 
electron at the Fermi energy to pass through the wire. In this case the Kondo temperature, Tk is unaffected by the 
finite wire length, L. On the other hand, when L <C £_r- the Kondo screening cloud must leak into the external leads. 
In fact, we might then think of the quantum wire plus dot as acting like a large quantum dot. (But note that, for 
simplicity, we haven't included an associated charging energy in the Hamiltonian.) 

Ignoring the Kondo interactions, the weak tunnelling between wire and leads, thw i leads to a broadening of the 
discrete spectrum of the wire. That is to say, the local density of states, ocurring in the Kondo interaction, has peaks 
at the energies given by the discrete spectrum of the wire, with spacing A = nvp/L and broadenned by an amount 
S oc t 2 LW . We assume that j « A. When L <C £rt, it was found that the conductance of the device depends very 
sensitively, at low temperatures, on the precise value of ew, which shifts the broadened energy levels of the wire, 
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FIG. 8: Quantum dot embedded in finite length quantum wire, eo and e^ control respectively the dot and wires gate voltage. 
Here we only consider the parity symmetric case whre all "L" and "_R" quantities are equal. 

relative to the Fermi energy, moving them on or off resonance. The renormalization of the Kondo coupling may be 
calculated straightforwardly in second order, using this density of states. Off resonance, there is little renormalization 
of the Kondo coupling down to exceedingly low energy scales. On resonance, there is again little renormalization until 
the effective band width is reduced to 6, the width of the resonance peak. At that stage the large density of states 
leads to a rapid growth of the effective Kondo coupling. In this situation it is convenient to define an effective Kondo 
temperature, Tk which depends on the length, L, of the wire and also on the gate voltage applied to the wire, ew- 
The non-resonant effective Kondo temperature, Tk when Tk 3> A, i.e. when the screening cloud fits inside 

the wire. However, T^^ - drops rapidly to zero as Tk is reduced below A. On the other hand, while the on-resonance 
Kondo temperature, T^, also tracks Tk when T 3> A, as Tk is reduced further, T^ levels off at a value of order 5, 
the width of the resonance controlled by t^w-, as illustrated in Fig. (|9]). Correspondly, the on-resonance conductance 
through the device, when S <C Tk <C A, stays very small down to temperatures of order S before shooting up towards 
the unitary limit, quite unlike the behaviour in the case A <C Tk, when the conductance approaches the unitary 
limit as the temperature is lowered below Tk- This is illustrated in Fig. (ITU|) . The dependence of the low T <C A 
conductance on the gate voltage, ew, applied to the wire, is also interesting. When the screening cloud doesn't fit 
inside the wire, L <C £,k, there are resonant peaks with spacing in ew given by the level spacing of the decoupled 
wires of length L, A. On the other hand, when the screening cloud fits inside the wire, £k "C L, the resonant peaks 
in low T conductance are separated by A/2, the level spacing for the doubled wire of length 2L. 



FIG. 9: The curves— represent Tk = f(T K ) in a log-log scale keeping the same values for A n ,S n for ew on resonance (plain 
curve which becomes almost flat at low T K <C 8 n ), and ew off resonance (dashed curve which drops sharply at low T K ). Both 
curves coincide at T K > A„ . 
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FIG. 10: On-resonance conductance 25 as a function of temperature in a symmetric device (assuming ew is on resonance) for 
both cases A n -C T% j (right blue curve) and A n 2> 5 n 2> T^ 2 e ft re d curve), from [2^|. The curves in plain style correspond 
to the perturbative calculations plus the Fermi liquid result for the first case only. We have schematically interpolated these 
curves (dotted lines) where neither the perturbative nor the Fermi liquid theory applies. 



D. The Kondo Box Model 



In 1999, Thimm, Kroha and van Delft introduced^ a model of a Kondo spin in a finite size system, which they 
referred to as "the Kondo box" model. It behaves quite differently than the models reviewed so far in this article, 
and here we focus on this difference. The experimental application that these authors discussed, a single magnetic 
impurity in a disordered ultra-small metallic grain, was rather different than those that we have discussed so far. 
They adopted a simplified Anderson impurity model in which: 

1. All energy levels of the unperturbed system (at Jk = 0) were assumed to be equally spaced in energy near the 
Fermi energy. 

2. All unperturbed states with energies near the Fermi energy were assumed to have equal hybridization with the 
impurity. 

The corresponding Anderson model is: 

H = ^£icj Q Cj a + eo^d^da + Uh d ^n di + v^^{c\ a d a + h.c). (129) 

ia a i,a 

where the e; are equally spaced. Taking the large U limit, and ignoring potential scattering, the corresponding Kondo 
model is: 

ia i,j 

where S is the spin of the magnetic impurity, S = Sad/2 and J = v 2 /U. The (uniform) energy level spacing, near 
the Fermi energy was estimated as: 

mkF V 

where V is the volume of the grain. Note that in the ID models considered in this section, the level spacing is 
proportional to 1/L, not 1/V. The behaviour of the Kondo box model follows immediately from the discussion of the 
ID models earlier in this section, if we make the corresponding replacement for A. As usual, finite size effects become 
important when Tfc ~ A, the bulk Kondo temperature is of order the level spacing. This implies^ that grains of order 
10 ran. in linear dimension are required to see finite size effects, considerably smaller than the estimate of 1 micron 
in the ID case. The behaviour of this model seems to contradict the naive Kondo screening cloud picture since there 
are essentially no finite size effects even when the linear dimension of the grain is much smaller than vf/Tr. 

In fact, the behaviour of this model depends strongly on the assumptions 1 and 2, above, which certainly don't 
always apply. To make this point it is instructive, as a thought experiment, to consider another model of an ultra-small 
grain which is far less realistic than even the above one. Imagine that 
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1. The grain is completely clean with no disorder 

2. The Kondo interaction is short range and spherically symmetric, the standard assumptions in the Kondo model 

3. The grain is perfectly spherical with the impurity atom exactly at the centre 

In this case, we may expand in spherical harmonics and project onto the s-wave, as usual. The s-wave model is 
essentially the one studied in sub-section III A| with the impurity spin at one end of a finite length ID wire. Now finite 
size effects set in when the length L is of order vf/Tk — £k as discussed above. The essential difference between 
the two models (or two sets of assumptions) is that in the second case most of the states of the metal have zero (or 
negligible) coupling to the impurity, whereas in the first case, all states are assumed to have equal coupling. In the 
second model only the s-wave states couple, a negligible fraction of the total set of states. This happens because the 
higher harmonic wave-functions all vanish at the origin, where the impurity is assumed to be located. Note also that 
the Kondo temperature will be much lower for this spherical model since the relevant density of states, appearing in 
the dimensionless Kondo coupling, scales with the volume in the first case but only as the length in the second. 

The behaviour of the second model is very strongly dependent on the assumptions, as can be seen by relaxing only 
the third of them. Suppose that we allow the surface of the grain to have a non-spherical shape. Then the usual 
wave-functions, involving products of spherical harmonics and spherical Bessel functions, are no longer eigenstates 
of the non-interacting model. If we insist on using this basis, then the impurity only interacts with the s-wave as 
usual. However, the reflections of the spherical harmonic wave-functions off the walls of the grain mix these wave- 
functions together. Actually, analogous phenomena occur in the ID models considered earlier in this section if we take 
more realistic models of the quantum wires with several channels ) 20,25 corresponding to transverse wave-functions. 
Consider, for example, the model in Sec. IIIIBI of a quantum dot at the end of a finite length quantum wire, which 
has several channels. The quantum dot hybridizes with some linear combination of channels in the wire. We can 
always adopt a basis such that one basis state is the one hybridizing with the quantum dot. However, these states 
are mixed by scattering off the other end of the wire (away from the quantum dot). A simple, intuitive picture, is 
that the screening cloud wave-function is now a linear superposition of the various channels, or the screening cloud 
folds back upon itself many times inside the wire, due to the reflections off the ends. If the wire has N channels, 
its effective length could be as long as LN. An important point is that embedding a Kondo spin in a finite system 
typically destroys one of the key simplifying features of the Kondo model - that the impurity spin only couples to one 
channel of unperturbed eigenstates (eg. the s-wave). The only case where we can be sure this assumption is still valid 
is when the finite system is an ideal single channel ID system. 

Depending on various assumptions about the grain shape, magnetic impurity location and disorder, one could 
imagine the Kondo temperature at which finite size effects set in ranging anywhere from order nvp/L to order 
7T 2 /m/EfV. In all cases, we expect that the condition for finite size effects to be important is that Tk be of order the 
finite size gap, suitably weighted by the hybridization of the corresponding wave- functions with the impurity. Only 
in single-channel situations does the simple Kondo screening cloud picture apply. 

Note that this complication makes observation of the screening cloud relatively easier in bulk situations where there 
are no boundaries to mix the hybridizing wave-function with the other channels. However, even there, it seems likely 
that scattering off other Kondo impurities might have a similar effect. 

IV. CONCLUSIONS 

In the standard Kondo model, describing a single magnetic impurity in an infinite host, the length scale = vp /Tk 
appears as a characteristic crossover length scale for many (perhaps essentially all) physical observables. Some intuition 
about the nature of the crossover is provided by imagining that, in the ground state, the impurity spin forms a spin 
singlet with a single conduction electron, whose wave-function is constructed out of states within energy Tk of the 
Fermi surface and hence is spread out over a distance of order One must be careful to only use this picture at 
low T <C Tk and even then it sould be used with caution. Naively, we might expect to see weak coupling behaviour, 
even at T = 0, at distances r <C S,k and strong coupling behaviour, outside the screening cloud, at r > (jf. While 
this seems to apply to the density oscillations and static impurity spin correlations, the first of these expectations 
turns out be be wrong, or at least to need some modifications, in the case of the Knight shift. The crossover effects 
at £k appear to be very hard to observe experimentally, precisely because they occur at such long distances. In a 
D-dimensional sample, they are suppressed by a factor of l/r D (or l/r ^ 1 ) making their observation more favourable 
in systems of reduced dimensionality. 

In a mesoscopic single-channel ID systems with a component of finite length, L, there is again typically a cross-over 
in physical quantities when £k is of order L. In all cases, weak coupling behaviour occurs when L -C £k so that the 
screening cloud doesn't "fit" inside the device and strong coupling behaviour occurs when £k -C L so that finite size 
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effects become unimportant. Except in ideai circumstances, mesoscopic systems do not have only a single channel. 
While this is unimportant for the Kondo effect in bulk, it can have important consequences for the observation of the 
Kondo screening cloud in a mesoscopic device. For an iV-channel system, the effective length "seen" by the Kondo 
screening cloud may be as large as LN . 

Throughout most of this paper we have ignored Coulomb interactions between the conduction electrons, only taking 
them into account at the magnetic impurity or quantum dot. While this may be justified at low T by Landau's Fermi 
liquid theory, for D=2 or 3, it would be necessary for the inelastic scattering length to exceed for our analyses to 
apply. In the D=l case, such interactions are important, even at T = 0. We analysed them only in a particularly 
simple case, the quantum dot end-coupled to a single channel quantum wire. While they somewhat reduced effects 
associated with the Kondo screening cloud, we concluded that the effects are still present and qualitatively unchanged. 

Another effect which can interfere with observation of the screening cloud is non-magnetic disorder. We only 
touched on this briefly in Sec. IIIID1 is our review of the "Kondo box" model. In addition a finite density of magnetic 
impurities further complicates the situation. 

Throughout essentially all of this review, we have considered only the Kondo model, not the more realistic Anderson 
impurity model. We expect similar crossover effects to also occur in that case but with another characteristic length 
scale entering, where the crossover from the free orbital to (unstable) local moment fixed point occurs. The Kondo 
crossover can only be observed if this other crossover occurs at a significantly shorter length than 

After many years of theoretical and theoretical investigations, the Kondo screening cloud, as predicted by the basic 
Kondo model, remains undetected experimentally. This is likely due to a combination of experimental difficulties, 
and limitations of the basic Kondo model. It is to be hoped that further experimental and theoretical progress will 
eventually bring this long search to a happy conclusion. 



Acknowledgements 

I would like to thank all my collaborators in this work. In chronological order: Erik S0rensen, Victor Barzykin, 
Pascal Simon, Rodrigo Pereria, Nicolas Laflorencie, Bert Halperin, Hubert Saleur and Lazlo Borda. I have also 
benefited from discussions on this topic with Jan von Delft, Junwu Gan, Nikolai Prokof'ev and Fred Zawadowski. 
This research was supported in part by NSERC and ClfAR. 



1 A.C. Hewson, The Kondo problem to heavy fermions, [Cambridge University Press, Cambridge (1993)]. 

2 I. Affleck and A.W.W. Ludwig, Nucl. Phys. B360, 641 (1991). 

3 K.G. Wilson, Rev. Mod. Phys. 47, 773 (1975). 

4 N. Andrei, Phys. Rev. Lett. 45, 379 (1980); N. Andrei, K. Furuya and J.H. Lowenstein, Rev. Mod. Phys. 55, 331 (1983). 

5 PB. Weigmann, JETP Lett. 31, 364 (1980). 

6 L. Borda, Phys. Rev. B75, 041307 (2007). 

7 I. Affleck, L. Borda and H. Saleur, Phys. Rev. B77, 180404 (2008). 

8 W.B. Thimm, J. Kroha and J. von Delft, Phys. Rev. Lett. 82, 2143 (2001). 

9 V. Barzkyin and I. Affleck, Phys. Rev. B61, 6170 (2000). 

10 V. Barzykin and I. Affleck, Phys. Rev. Lett. 76, 4959 (1996). 

11 V. Barzykin and I. Affleck, Phys. Rev. B57. 432 (1998). 

12 J. Gan, J. Phys.:Cond. Mat. 6, 4547 (1994). 

13 E. S0rensen and I. Affleck, Phys. Rev. B53, 9153 (1995). 

14 J.P Boyce and CP Slichter, Phys. Rev. Lett. 32, 61 (1974); Phys. Rev. B13, 379 (1976). 

15 G. Bergman, Phys. Rev. B77, 104401 (2008). 

16 F. Mezei and G. Griiner, Phys. Rev. Lett. 29, 1465 (1972). 

17 F. Lesage and H. Saleur, J. Phys. A30, L457 (1997). 

18 H. Ishii, J. Low Temp. Phys 32, 457 (1978). 

19 I. Affleck and P. Simon, Phys. Rev. Lett. 86, 2854 (2001). 

20 P. Simon and I. Affleck, Phys. Rev. B64, 085308, 2001. 

21 E.S. S0rensen and I. Affleck, Phys. Rev. Lett. 94, 086601 (2005). 

22 H. Frahm and A. A. Zvyagin, J. Phys. 9, 9939 (1997). 

23 R.G Pereira, N. Laflorencie, I. Affleck and B.I. Halperin, Phys. Rev. B77, 125327 (2008). 

24 PS. Cornaglia and C.A. Balseiro, Phys. Rev. B66, 115303 (2002); ibid 174404 (2002); Phys. Rev. Lett. 90, 216801 (2003). 

25 P. Simon and I. Affleck, Phys. Rev. Lett. 89, 206602 (2002); Phys. Rev. B68, 115304 (2003). 



